August 20 



We present an improved direct measurement of the parity-violation parameter At, in 
the Z boson-b quark coupling using a self-calibrating track-charge technique applied to 
a sample enriched in Z — > bb events via the topological reconstruction of the B hadron 
mass. Manipulation of the SLC electron-beam polarization permits the measurement of A b 
to be made independently of other Z-pole coupling parameters. From the 1996-98 sample of 
400,000 hadronic Z decays, produced with an average beam polarization of 73.4%, we find 
A b = 0.906 ± 0.022(stat.) ± 0.023(syst.). 
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Measurements of b quark production asymmetries at the Z° pole determine the extent 
of parity violation in the Zbb coupling. At Born level, the differential cross section for the 
process e + e~ — > Z° — > bb can be expressed as a function of the polar angle 9 of the b quark 
relative to the electron beam direction, 

a (cos 8) = dab /d cos 6 oc (1 - A e P e )(l + cos 2 9) + 2A b (A e - P e ) cos (9, (I) 

where P e is the longitudinal polarization of the electron beam (P e > for predominantly 
right-handed polarized beam). The parameter Af = 2vfdf/(vf + aj), where Vf (a/) is the 
vector (axial vector) coupling of the fermion / to the Z° boson, with / = e or b, expresses 
the extent of parity violation in the Zff coupling. 

From the conventional forward-backward asymmetries formed with an unpolarized elec- 
tron beam (P e = 0), such as that used by the CERN Large Electron- Positron Collider 
(LEP) experiments, only the product A e Ab of parity-violation parameters can be measured 
[1]. With a longitudinally polarized electron beam, however, it is possible to measure Aj, 
independently of A e by fitting simultaneously to the differential cross sections of Eq. (1) 
formed separately for predominantly left- and right-handed beam. The resulting direct mea- 
surement of Ab is largely independent of propagator effects that modify the effective weak 
mixing angle, and thus is complementary to other electroweak asymmetry measurements 
performed at the Z° pole. 

In this Letter, we present a measurement of A b based on the use of an inclusive vertex 
mass tag (improved relative to that of previous publications due to the use of an upgraded 
vertex detector) to select Z — > bb events, and the net momentum-weighted track charge [2] 
to identify the charge of the underlying quark. This result, incorporating data collected 
during the 1996-98 runs of the Stanford Linear Collider (SLC), is over twice as precise as 
that of our previous publication [4], which was based on data from 1993-95. 

The operation of the SLC with a polarized electron beam has been described elsewhere [5]. 
During the 1996-98 run, the SLC Large Detector (SLD) [6] recorded an integrated luminosity 
of 14.0 pb -1 , at a mean center-of-mass energy of 91.24 GeV, and with a luminosity-weighted 
mean electron-beam polarization of \P e \ = 0.734 ± 0.004 [3]. The 1996-98 run of the SLD 
detector incorporated the upgraded VXD3 CCD pixel vertex detector [7], which featured a 
greater coverage in cos 9, as well as a larger outer radius and substantially less material per 
layer, than that of the VXD2 vertex detector [8] in place from 1993-95. 

The SLD measures charged particle tracks with the Central Drift Chamber (CDC), which 
is immersed in a uniform axial magnetic field of 0.6T. The VXD3 vertex detector provides 
an accurate measure of particle trajectories close to the beam axis. For the 1996-98 data, 
the combined r<p (rz) impact parameter resolution of the CDC and VXD3 is 7.7 (9.6) fim 
at high momentum, and 34 (34) /iin at pj_Vsin9 = 1 GeV/c, where p± is the momentum 
transverse to the beam direction, and r (z) is the coordinate perpendicular (parallel) to the 
beam axis. The combined momentum resolution in the plane perpendicular to the beam 
axis is 5pj_/p± = \J (.01) 2 + (.0026 p±/GeV/c) 2 . The thrust axis is reconstructed using the 
Liquid Argon Calorimeter, which covers the angular range | cos#| < 0.98. 

The details of the analysis procedure are similar to those of the 1993-95 sample analysis. 
Events are classified as hadronic Z° decays if they: (1) contain at least seven well-measured 
tracks (as described in Ref. [6]); (2) exhibit a visible charged energy of at least 20 GeV; 
and (3) have a thrust axis polar angle satisfying | cos 9 t hrust\ < 0.7. The resulting hadronic 
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sample from the 1996-98 data consists of 245,048 events with a non-hadronic background 
estimated to be < 0.1%. 

We select against multi-jet events in order to reduce the dependence of the measured 
value of A b on the effects of gluon radiation and inter-hemisphere correlation. Events are 
discarded if they are found to have four or more jets by the JADE jet-finding algorithm with 
Ucut = 0.02 [9], using reconstructed charged tracks as input. In addition, any event found to 
have three or more jets with y cut = 0.1 is discarded. 

To increase the Z° — > bb content of the sample, a tagging procedure based on the invariant 
mass of 3-dimensional topologically reconstructed secondary decay vertices is applied [10]. 
The mass of the reconstructed vertex is corrected for missing transverse momentum relative 
to the reconstructed B hadron flight direction in order to partially account for neutral 
particles. The requirement that the event contain at least one secondary vertex with mass 
greater than 2 GeV/c 2 results in a sample of 36,936 candidate Z° — > bb decays. The purity 
(97%) and efficiency (77%) of this sample are calculated from the data by comparing the 
rates for finding a high mass vertex in either a single or both hemispheres, where the two 
hemispheres are defined relative to the plane perpendicular to the thrust axis. This procedure 
assumes a-priori knowledge of the small udsc tagging efficiency, as well as the size of inter- 
hemisphere correlations, both of which are taken from Monte Carlo (MC) simulation. This 
procedure also assumes knowledge of the Z — > cc and Z — > bb branching fractions, which 
are assigned their Standard Model values of 0.172 and 0.216, respectively. 

We construct a signed thrust axis T, which provides an estimate of the direction of the 
negatively-charged b quark, as follows. Using all track-charge quality tracks, as defined in 
Ref. [11], we form the track-direction-signed (Q) and unsigned (Q+) momentum- weighted 
track-charge sums 

Q = ~ E q r sMp r f)\(p r f)r, (2) 

tracks 

Q + = E qM-tT, (3) 

tracks 

where qj and pj are the charge and momentum of track j, respectively. T is chosen as 
the unit vector parallel to the thrust axis that renders Q > 0. We use k = 0.5 to max- 
imize the analyzing power of the track charge algorithm for Z° — > bb events, resulting in 
a correct-assignment probability of 70%. Fig. 1 shows the T z = cosOthrust distribution of 
the 6-enriched sample separately for left- and right-handed electron beams. Clear forward- 
backward asymmetries are observed, with respective signs as expected from the cross-section 
formula in Eq. 1. 

The value of A b is extracted via a maximum likelihood fit to the differential cross section 
(see Eq. 1) 

p\A b ) = (1 - A e Pl){l + (TIT) + 2(A e - Pl)T z [A b fl(2pl - 1)(1 - A* QCA6 ) + 
Acf c Wc - 1)(1 - Xqcd,c) + A bck9 (l -f b - f c )(2p\ ckg - 1)], (4) 

where P e * is the signed polarization of the electron beam for event i, fL c \ the probability that 
the event is a Z° — > bb(cc) decay (parameterized as a function of the secondary vertex mass), 
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and A l QCD,b,c are final-state QCD corrections, to be discussed below. Ab c k g is the estimated 
asymmetry of residual uu, dd, and ss final states. The parameters p are estimates of the 
probability that the sign of Q accurately reflects the charge of the respective underlying 
quark, and are functions of \Q\, as well as the secondary vertex mass and \T Z \. 

As in our previous publication [4], we measure pb directly from the data [12]. Defining 
Qb (Qi) to be the track-direction-unsigned momentum-weighted track charge sum for the 
thrust hemisphere containing the b (b) quark, the quantities 

Qsurn = Qb + Qb i Qdif = Qb ~ Qb i (5) 

may be related to the experimental observables defined in Eqs. 2 and 3 respectively: \Qdif\ — 
\Q\ and Q sum = Q+ ■ Our MC simulation indicates that the Qb and Qi distributions are 
approximately Gaussian. In this limit [12], 



with 



Hif _ *y/(\Q*f\ 2 )-<r; 
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dif 



ab = _^L = «/ , (7) 

a dif a dif 

where g^/ an d Cdif are the mean and width, respectively, of the Gaussian Qdif distribution. 
The parameter a b , whose magnitude depends upon the separation between the b and b track- 
sum distributions via the observable (\Qdif\ 2 ), provides a measure of the analyzing power 
of the 6-quark direction estimator T. Figure 2 compares the distributions of the observable 
combinations \Qdif\ and Q + between data and MC. 

In the absence of a correlation between Qb and Qj;, a dif = cr sum , where a sum is the ob- 
served width of the Q + distribution. Thus ab can be derived from experimental observables. 
In the presence of a correlation, adif — (1 + A)a snm , where A characterizes the strength of 
the correlation, which can be determined from the MC simulation. For JETSET 7.4 [13] 
with parton shower evolution, string fragmentation, and full detector simulation, A is found 
to be 0.040. The effects of light flavor contamination are taken into account by adjusting 
the observed widths o 2 sum and (\Qdif\ 2 ), using the magnitude and width of the light-flavor 
and cc contributions estimated from the MC. This correction increases the value of a& by 
2% to 0.2944 ± 0.0078, bringing it into good agreement with the value of 0.2949 ± 0.0007 
extracted from the Z — > bb simulation. 

Final-state gluon radiation reduces the observed asymmetry from its Born-level value. 
This effect is incorporated in our analysis by applying a correction Aqcd(\ cos 6*1) to the 
likelihood function (Eq. 4). Calculation of the quantity Aq C d has been performed by several 
groups [14]. 

For an unbiased sample of bb events, correcting for final-state gluon radiation increases 
the measured asymmetry by ~ 3%. However, QCD radiative effects are mitigated by the 
use of the thrust axis to estimate the 6-quark direction, the Z° — > bb enrichment algorithm, 
the self-calibration procedure, and the cut on the number of jets. A MC simulation of 
the analysis chain indicates that these effects can be represented by a cos ^-independent 
suppression factor, xqcd = 0.074, such that Aq C d = x qcd^qcd- 
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Effects due to gluon splitting to bb and cc have been estimated by rescaling the JETSET 
simulation production of such quark pairs to current world-average gluon splitting measure- 
ments [15], leading to a correction of +0.3% on the value of A b . Additional radiative effects, 
such as those due to initial-state radiation and 7/Z interference, lead to a further correction 
of —0.2% to the measured value of A b . 

While, as described above, the overall tagging efficiency is derived from data, the depen- 
dence of the 6-tagging efficiency upon the secondary vertex mass must be estimated from 
the MC simulation, as must be the charm correct-signing probability p c . The value of A c is 
set to its Standard Model value of 0.67, with an uncertainty commensurate with that of [16]. 
The value of A bckg is set to zero, with an uncertainty corresponding to the full physical range 
l^bcfcj < 1- The resulting value of A b extracted from the fit is A b = 0.907 ± 0.022 (stat). 
This result is found to be insensitive to the value of the ft-tag mass cut, and the value of 
weighting exponent k used in the definition (2) and (3) of the momentum-weighted track 
charge sum. 

We have investigated a number of systematic effects which can change the measured 
value of A b ; these are summarized in Table I. The uncertainty in a b due to the statistical 
uncertainties in (\Qdif\ 2 ) and af um corresponds to a 1.6% uncertainty in A b . The uncertainty 
in the hemisphere correlation parameter A is estimated by varying fragmentation parameters 
within JETSET 7.4, and by comparison with the HERWIG 5.7 [17] fragmentation model. 
The resulting uncertainty in A b is 1.4%. The sensitivity of the result to the shape of the 
underlying Q b distribution is tested by generating various triangular distributions as well 
as double Gaussian distributions with offset means. The test distributions are constrained 
to yield a Q sum distribution consistent with data, and the total uncertainty is found to be 
0.8%. In addition, while the mean value of the self-calibration parameter a b is constrained 
by the data, it has a cos 6 dependence due to the fall-off of the tracking efficiency at high 
I cos#| which must be estimated using the simulation, leading to a 0.4% uncertainty in A b . 

The extracted value of A b is sensitive to our estimate of the Z° — > cc background, which 
tends to reduce the observed asymmetry due to the positive charge of the underlying c quark. 
The uncertainty in the purity estimate of 96.9 ± 0.3% is dominated by the uncertainties in 
the charm tagging efficiency (e c = 0.0218 ± 0.0004) and the statistical uncertainty of the 
bottom tagging efficiency determined from data, leading to a 0.5% uncertainty in A b . An 
outline of the charmed quark efficiency uncertainty determination can be found in Ref. [18]; 
the uncertainty is dominated by empirical constraints on charmed hadron production rates 
and on K° production in the decay of charmed mesons. Uncertainties in the measured values 
of R b and R c contribute, through the tag purity, to uncertainties in A b of 0.1% and 0.0%, 
respectively. 

Agreement between the data and MC simulation charged track multiplicity distributions 
is obtained only after the inclusion of additional ad-hoc tracking inefficiency. This random 
inefficiency was parameterized as a function of total track momentum, and averages 0.4 
charged tracks per event, leading to an overall change of +1.3% in A b . As a check, we 
employ an alternative approach, matching the efficiency of the linking of the independent 
CDC and VXD3 track segments between data and MC simulation. This yields a change of 
+0.5% in A b ; we take the difference of 0.8% as an estimate of the systematic error on the 
modeling of the tracking efficiency. Combining all systematic uncertainties in quadrature 
yields a total relative systematic uncertainty of 2.6%. 
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The extracted value of A b depends on a number of model parameters, as follows. Increases 
by 0.01 in the values of A c , R b , R c , and the per-event rate of bb production via gluon splitting, 
lead to changes in A b of +0.0002, -0.0055, +0.0002, and +0.0110, respectively. 

In conclusion, we have exploited the highly polarized SLC electron beam and precise 
vertexing capabilities of the SLD detector to perform a direct measurement of A b = 0.906 ± 
0.022(stat) ±0.023(syst), from the 1996-98 SLD data sample. Combined with our previously 
published result [4] based on the 1993-95 data sample, we find 

A b = 0.907 ± 0.020(stat) ± 0.024(syst), (8) 

for the full 1993-98 data sample. This result is in good agreement with the Standard Model 
prediction of 0.935, and represents an improvement of over a factor of two in the precision 
of the determination of A b via the use of momentum-weighted track charge. 

We thank the staff of the SLAC accelerator department for their outstanding efforts 
on our behalf. This work was supported by the U.S. Department of Energy and National 
Science Foundation, the UK Particle Physics and Astronomy Research Council, the Istituto 
Nazionale di Fisica Nucleare of Italy and the Japan-US Cooperative Research Project on 
High Energy Physics. 
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FIG. 1. Polar angle distributions for track-charge-signed Z — > 66 candidates, separately for left- 
and right-handed electron beam. The shaded histogram represents the contribution from non-66 
background, estimated as described in the text. The analysis employs a cut of | cos#| < 0.7. 
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TABLES 



Error Source 


Variation 


<r A I a 

5A b /A b 


Self- Calibration 






a b statistics 


±1(7 


1.6% 


Xb Correlation 


JETSET, HERWIG 


1.4% 


P(Qb) shape 


Different shapes 


0.8% 


cos 6 shape of a b 


MC Shape vs Flat 


0.4% 


Light Flavor 


50% of correction 


0.2% 


Analysis 






Tag Composition 


Procedure from [18] 


0.5% 


Detector Modeling 


Compare tracking 


0.8% 




eff. corrections 




Beam Polarization 


±0.5% 


0.5% 


QCD 


Full correction 


0.3% 


Gluon Splitting 


Full Correction 


0.1% 


A c 


0.67 ±0.04 


0.1% 


Abckg 


0±0.50 


0.2% 


Total 




2.6% 



TABLE I. Relative systematic errors on the measurement of A b . 
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